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On a modified Hulthén potential 


By S. G. M. Gusravi 


Introduction 


Meson theory has long ago predicted highly singular potentials of the form 
e/(ur), e “/(ur)? and even e “’/(wr)? for the two-nucleon interaction. Because 
of analytical difficulties these potentials have been treated almost exclusively with 
numerical methods, especially when they are used with a “‘cut off’? which removes 
the singularity, as is usual and in some cases necessary. Being forced to numerical 
solutions, however, one loses much of the general survey of the problem and therefore 
it may be of some value to construct modified potentials, similar to those above in 
form and physical properties, which can be treated analytically. For the Yukawa 
potential we already have the well-known Hulthén potential e~”/(1—e~“") [1], and 
here we shall discuss a natural modification of this potential to the forme “"/(1—e “")?, 
which is as easy to handle analytically as the original one and which may be used 
as an approximation to the meson potential e“”/(wr)?. At least the singular properties 
at the origin must be the same for the approximative and the original potential and 
for large values of ~r where the approximation is not very good both potentials are 
so small in value as to be neglected. An analogous approximation for the potential 
e-“"/(ur)? is of course also conceivable but in this case it seems as if the approximative 
potential does not have any special advantage over the original one. 


Solution of the wave equation 


In the S state which is the only one of interest here the radial wave equation for 
the deuteron system may be written in common notation 


eat + GB VOn)] RO)=0, a) 

where M is the nucleon mass, # the energy and V(r) the (central) potential. Here 
ats 

Be Fo ms’ (2) 


where V, is a positive constant. After the usual transformation R(r) =w(r)/r we 


get from (1): 7 
E ao | u(r) = a V (r) u(r), (3) 


31 437 


Ss. G. M. GUSTAVI, On a modified Hulthén potential 


where k? = ME/h. In the case of an ordinary Hulthén potential this equation can | 
be transformed into a hypergeometric differential equation by the transformation 


enn (eo) = eM =z 


(see [1] and [2]), and this is suited also for our problem. After applying it to equation | 
(3) with V(r) taken as the modified Hulthén potential (2) we get 


a? 2ik\l dd % 
— F (z)=0, 4)) 
gat (1+ lt rata (2) (9 


where v, = (M/h?)(Vo/u2). This is not a hypergeometric equation, but it is, however, . 
a second-order differential equation with three regular singular points and can thus: 
be transformed into a hypergeometric equation [3]. The suitable transformation is : 


F(z) = (1—z)' f (2), (5) 
where ¢ is one of the roots of the indical equation for the singular point z = 1 (r = 0): : 
t(é—1)+- 95 = 0: (6) | 


We shall in the following denote these two roots by 


t=44+V4—-%, t,=4—-V}—%, (7) 


while ¢ without index may represent any of them. The transformation (5) will bring 
(4) into the form 


ete— yy" (e)+ |(1+2¢+2"*) 2—(1+224}] f (2)+ Je +20] f)=0. (8) 
le ia le 


This is a hypergeometric equation, the general form of which may be written [4]: 


2(z—1) 9" (z) + [(1 +a + 6)z—c] g’(z) + ab g(z) =0. (9) 
Comparison between (8) and (9) gives 
2; ; 
a=t; p eepptsl Se ging oth (10) 
ee & 


The solutions of (9) are most suitably given in terms of series solutions about 
some or other of the three singular points of the equation, and all these series solu- 
tions can be expressed with help of the hypergeometric series [4]: 


F(a, b,c; Hela ea : 
c 


|2e(c+1) 
__ Te) ST (a+n)l(b+n) , 
Pa) P Oi Sea Li(eaan) Pilea & (ul) 


Weare primarily interested in two kinds of solutions of (3) (and correspondingly of 
(8)): for the discussion of bound state problems we need the negative energy (ik =x, 
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x real and positive) solutions that are regular at infinity and for low energy scattering 
problems we need the zero energy (k = 0) solutions. 
We begin with the first case. The general series solution of (9) about z =0 (7 = 09) 
is [4]: 
g(z) =A F (a,b,c; z) + Bz’ ° F(a—c + 1,b-—c +1, 2—c; 2), 


which gives according to (10): 


9), 2x 


2; 2% / Ox Qn | 
j= AP (042% 1+="12) + Be * w(t oe Ay 2) (12) 
lu lu | 


From this it is easy to see that the solution of the wave equation (3) for ik = x(x real 
and positive) that is regular at infinity (for 7 =°co or z=0) must be of the form 


24 ee 
> © 


wie —e mye (e+ , l+ (13) 
be [ 


This solution at first sight seems to depend on the choice of ¢ in an essential way 
but in fact it does not. In order to realize that and to obtain a form of the solution 
which is sometimes more useful than (14) we employ the analytic continuation for- 
mula [4] 

I(c) ' (c—a—6b) 
I'(c—a)T' (c—}) 


['(c) (a+ 6—c) 
T'(a) 1 (0) 


Ff (a, b,c; z)= F(a,6,a+b—c+1; 1—2)+ 


(l—z)° °°? F(c—a, c—b, c—a—b+1; 1-2), (14) 


a 


which inside the common region of convergence connects the ordinary hypergeo- 
metric series with the series solutions of equation (9) about z = 1. (13) and (14) give 


byt, 2% 

TEN WAS Pee 

2 2 lu 

) 2% 
a Bea eek cw R (ty a i, oad (15) 
2x 

it (1.422) 
Pat (aor 


where we have taken ¢ = ¢, in (13) and used the relations 
t, +t, =1 and [(z+1)=2zT (2). 


Apart from a minus sign the choice ¢ =f, in (13) will obviously give the same result. 
It may be remarked that (15) shows u(x; r) to be real also in the case of complex ¢ 
values (v,) > 4), because ft, is then the complex conjugate of f,. 

For wr<1 the hypergeometric series in (15) are approximately constant while 
(1 —e~"")*' = (ur)" and (1 — e '")** ~ (wr). Now from (7) if vy > 0 we get 0 < Re(t)<1. 
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Consequently w(x; r) always goes to zero with r, while the original radial function | 
R(x; 1r) =u(x;7)/r is always singular at the origin. If the ¢ values are real, that is \f 
if vp <+ the general form of the wave function is rather simple; it is always positive | 
and the second derivative changes sign at one place only. In the case of complex ¢ | 
values (v, > }), on the other hand, the wave function oscillates with a frequency 
which goes to infinity as r>0, because of factors of the form cos(v log ur), where y | 
is the imaginary part of t, (or f,). The distinction between real and complex values is | 
a very important one and we shall in the following denote the potential as wnder- 
critical or overcritical if vy< } or >} respectively. 

We now turn to the zero energy case (k =0). The general series solution of (9) 
about the singular point z = 1 (which corresponds to r =0) is: 


g(z) =AF(a,b,a+b—c+1;1—z)+ B(1—2z)"*’ F(c—b,c—a,e—a—b +1; 
1 =z); 
For k = 0 this gives two independent solutions of the equation (3) which may be written 


P(r) = (Loe) F(t, hy 2445 1 ane 


u 
ug (r)=(1—e )* (fo; te, Stale © )e 7) 
For small r values the behaviour of the functions above is quite similar to that of 
u(x; r) and especially we can distinguish between an undercritical and an overcritical 
case in the same way. In the overcritical case the functions (17) are not real, however, 
and it is then perhaps more natural to take the real and imaginary parts of one of 
them as independent solutions. 
The asymptotic behaviour of the zero energy functions for large 7 values needs 
some special attention. Let us define a, according to 


E(t. ¥5 265 2) 


Te) 3 Pete) *_ Ss 
M(@) oT (ttn) [2 wo” 


A straightforward estimation with help of the Stirling formula for the gamma 
function will yield the following formula for the asymptotic behaviour of an: 


Lej ARs 1 
oo Tay (ana as) 8 


Let us denote (1/n) [I (2 t)/I()] by 6, and consider the expression 


(lest) + & (anda) (1 aD +(l-e")' Sd (-e"y" (19) 
n= n=1 


which is identical with uq)(r) if we choose t = ¢; (i = 1, 2). Because 


Sil 
2 wine oa ~ lee (le ape 
we have lim Pin r—(l— —uryt ~ b —pryn | 
poe [Tl han he estate ]-o. (20) 
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The first term in (19), on the other hand, tends to a constant as r—>oo, because, 


according to (18), 
(28) (# 1 
cn = ba On = Tara ([+0(-)}, (n=1) (21) 


nN nv 


so that > c, is convergent and we may use Abel’s theorem for power series [5] so as to 
obtain 


lim Jaen (1+ S (ay br) 1 e-"7))| 
n=1 


Tow 


. [(2t) S /T?(t+n) 1 1 wi 
n=1 I’ (t) > (rare n) |n = (22) 


where by C(t) we have indicated that the constant depends on tf. It is easily seen that 
C(1) =0 and jim C(t) = 1, while between these values C(t) is a monotonously decreas- 


ing function of ¢, which must be evaluated numerically by summing the series (22). 
Putting together the formulae (19), (20) and (22) we get at last the asymptotic 
expression 


UD (r) ~ (t,) + = ih) ur. (23) 


Before we leave this discussion of the solutions of the differential equation (3) we 
should perhaps point out that in the case t, = t, = } some of the formulae above must 
be modified. In this case the two series about z = 1(7 =0) which are represented in 
(15) and (17) are identical and a second independent solution must be found in the 
form of a series which contains a logarithmic term and gives a wave function w(7r) 
with a behaviour of the form )/r log r at the origin. Consequently the formula (15) 
is not valid for t, = t, = 4 and one of the functions (17) must be replaced by the above- 
mentioned other series. 


Physical consequences 


The singularity of the modified Hulthén potential at the origin is in some respect a 
limiting case of what may be allowed. As we have seen all solutions of the radial 


wave equation (3) goes to zero at the origin either as r’ or r** (or Vr log r for the 
special case t, =t, = 4) while the radial function R(r) =u(r)/r is always singular at 
the origin. This last fact is not a reason for us to reject the potential (in the Dirac 
theory of the hydrogen atom we have a similar case), but on the other hand it is 
certainly not possible to accept both 7” and r* as allowed forms of the wave function 
for small 7, because that would destroy the uniqueness of the wave function for positive 
energies and give us a continuous spectrum for negative energies. As has sometimes 
been observed in the literature (see for instance [6], where scattering by an inverse 
cube law force is briefly discussed) the only way of using the potential in physical 
problems is to accept one of the ¢ values and reject the other one. But this would 
perhaps not be more than a mathematical device if we could not find any physical 
support for the procedure. 

Now, in the overcritical case (vj > 4) it is certainly impossible to make any physical 
distinction between the two exponential values and thus we must conclude that an 
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overcritical modified Hulthén potential can be used with cut off only. We may, 
however, treat the undercritical case independently of the overcritical one, although 
this is perhaps an ambiguous point. When the exponentials are real the “degree of 
singularity” in the wave function R(r) is not the same for both values and it may 
seem natural to accept the largest of them (corresponding to ¢,) and reject the other 
one. From a physical standpoint this choice is also supported for instance by a criter- 
ion used earlier in the relativistic treatment of the hydrogen atom referred to above 
(see for instance [7]) and which states that if a potential gives rise to a singularity in the 
wave function, this singularity must vanish together with the potential if it is to be 
accepted; this is only the case for the exponential ¢ =¢,, which tends to | as v) ap- 
proaches zero. 

To summarize the situation, we may state that complex exponentials and the 
exponential ¢, in the real case are unpermissible for physical reasons, while the real 
exponential t, perhaps may be accepted; in any case it is natural to examine the 
consequences of an acceptance. 

Let us first discuss the possibility of bound states for an undercritical potential. 
If the singularity arising from (1 — e-“”)* is not allowed, the factor I (2#,)/T' (4) T(t, + 
+ 2x/u) in the solution (15) must vanish for a bound state. But in the undercritical 
case all arguments in these gamma functions are real and positive, thus also the func- 
tions themselves (for t, =t, = 4 we get instead a non-vanishing part of the function, 


which behaves as )r log r at the origin). Consequently there are no bound states and, 
as it ought to be, this conclusion is still valid if we introduce a cut off in the potential, 
because the wave function u(x; 7) is always positive outside the origin if vy < +. The 
nonexistence of a bound state for an undercritical potential is also consistent with 
the fact that the scattering length is negative in this case as is evident from the asymp- 
totic form (23) of the zero energy wave function. 

We then turn to the low energy scattering problem. The allowed zero energy solu- 
tion is (17) with ¢ = ¢, and the asymptotic form of this, normalized as in the effective 
range theory is obtained from (23): 


k (t,) . k(t,) 7 P 
Tey or lbiouyt pe (24) 


where r; is the range of the potential and k(t) a notation for I'(2t)/T2(¢). If we want 
to adjust the undercritical potential to known singlet scattering data, (24) is to be 
compared with the form 1 —7r/a;, where we may take the experimentally given 
scattering length a; = — 23.69 x 10-13 cm [8]. For a fixed value of ¢,, this will give 
us an expression for 75: 


ie? (25) 
1 


and it is easy to see that this is a monotonously increasing function of t,. Now the 
smallest possible value of ¢, is 4 and the corresponding value of ry, computed from 
(25), is rs; = 8.51x 10-18 cm. But already a very rough estimation shows that this — 
value is much too large to give an effective range of the correct order of size ( 
= 2.44 0.5 x 10-38 cm according to [8]). 

If we introduce a cut off on the other hand, it is not difficult to show that in the 
undercritical case the potential is still too weak to be adjusted even to the singlet 


Tos = 
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scattering data, and in fact it is weaker than the corresponding potential without 
cut off as it ought to be. 

The somewhat disappointing result of our investigation of the undercritical 
modified Hulthén potential is consequently that this potential alone is too weak to 
be of any interest for the two-nucleon problem, and this result is certainly also valid 
for the original meson potential e “"/(ur)?. The investigation has shown, however, 
that the choice of t, as the only allowed exponential, gives physical results that are 
quite analogous to those for a more regular potential. 

An overcritical potential, which as we have seen, must be used with cut off, is 
of cause possible to apply both to the triplet and the singlet case, but unfortunately 
the wave functions are much more difficult to handle analytically if a cut off is intro- 
duced and we shall not here discuss this point any further. 

There still remains the possibility, however, to combine an ordinary and an 
(undercritical) modified Hulthén potential, that is to use a potential of the form 


— er —pr 


é 
(Ucn) 


é 
(l—e"? 


Var y= SV, V, 


where V, and JV, are constants. In this case also, the wave equation is soluble in 
terms of hypergeometric functions and the singular properties at the origin are quite 
the same as for the modified potential alone. Up till now, no detailed calculations 
have been carried out on the basis of this combined potential, however. 
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